OPEN 



ORIGINAL ARTICLE 



The ISME Journal (2013) 7, 1651-1660 

© 2013 International Society for Microbial Ecology All rights reserved 1751-7362/13 
www.nature.com/ismej 



Microbial biodiversity in glacier-fed streams 



Linda Wilhelm 1 , Gabriel A Singer 1,2 , Christina Fasching 1 , Tom J Battin 1 ' 2 and 
Katharina Besemer 1 ' 2 

^Department of Limnology, University of Vienna, Vienna, Austria and 2 Wasser Cluster Lunz, 
Lunz am See, Austria 



While glaciers become increasingly recognised as a habitat for diverse and active microbial 
communities, effects of their climate change-induced retreat on the microbial ecology of glacier-fed 
streams remain elusive. Understanding the effect of climate change on microorganisms in these 
ecosystems is crucial given that microbial biofilms control numerous stream ecosystem processes 
with potential implications for downstream biodiversity and biogeochemistry. Here, using a space- 
for-time substitution approach across 26 Alpine glaciers, we show how microbial community 
composition and diversity, based on 454-pyrosequencing of the 16S rRNA gene, in biofilms of 
glacier-fed streams may change as glaciers recede. Variations in streamwater geochemistry 
correlated with biofilm community composition, even at the phylum level. The most dominant phyla 
detected in glacial habitats were Proteobacteria, Bacteroidetes, Actinobacteria and Cyanobacteria/ 
chloroplasts. Microorganisms from ice had the lowest a diversity and contributed marginally to 
biofilm and streamwater community composition. Rather, streamwater apparently collected 
microorganisms from various glacial and non-glacial sources forming the upstream metacommu- 
nity, thereby achieving the highest at diversity. Biofilms in the glacier-fed streams had intermediate a 
diversity and species sorting by local environmental conditions likely shaped their community 
composition, a diversity of streamwater and biofilm communities decreased with elevation, possibly 
reflecting less diverse sources of microorganisms upstream in the catchment. In contrast, |! 
diversity of biofilms decreased with increasing streamwater temperature, suggesting that glacier 
retreat may contribute to the homogenisation of microbial communities among glacier-fed streams. 
The ISME Journal (2013) 7, 1651-1660; doi:10.1038/ismej.2013.44; published online 14 March 2013 
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Introduction 

Glacial ecosystems, at the interface between 
the cryosphere, hydrosphere, pedosphere and the 
biosphere, are particularly prone to impacts of 
climate change (Brown et ah, 2007a; Milner et ah, 
2009; Finn et ah, 2010; Woodward et al„ 2010; 
Jacobsen et ah, 2012). Glacier-fed streams constitute 
a prominent geomorphological and ecological com- 
ponent of the glacier foreland, and integrate 
upstream catchment processes. Changes in runoff 
following glacier retreat shift the relative contribu- 
tions of icemelt, snowmelt and groundwater to 
stream discharge (Brown et ah, 2007a; Milner 
et ah, 2009), thereby influencing streamwater tem- 
perature, electrical conductivity and channel stabi- 
lity in glacier-fed streams (Hannah et ah, 2007; 
Brown et ah, 2007a; Milner et al. 2009). The close 
hydroecological coupling between glaciers and 
glacier-fed streams makes biodiversity in these 
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streams particularly susceptible to glacier retreat 
(Brown et al., 2007a; Finn et al, 2010; Jacobsen 
et al., 2012). For instance, reduced glacial runoff 
accompanying glacier retreat facilitates the displa- 
cement of cold-adapted invertebrate species by 
range expansion of less cryophilic species towards 
higher elevations (Brown et al., 2007a; Milner et al., 
2009; Finn et al, 2010). Ultimately, such shifts in 
species distribution may threaten the biodiversity in 
isolated high-elevation ecosystems with yet 
unknown implications for downstream biodiversity 
(Hughes et al., 2009; Finn et al., 2011). 

Ice ecosystems harbour complex microbial com- 
munities (Simon et al., 2009; Anesio and Laybourn- 
Parry, 2012), which have been proposed as sentinels 
of climate change and the general attrition of the 
cryosphere (Vincent, 2010). The most abundant 
bacterial phyla detected in glacier ice are Proteo- 
bacteria, Actinobacteria, Bacteroidetes, Firmicutes 
as well as photosynthetic Cyanobacteria (Simon 
et al., 2009; Xiang et al., 2009; Anesio and Laybourn- 
Parry, 2012; Edwards et al, 2011). While the 
microbial diversity and activity in glacier ice was 
the focus of several recent studies (Hodson et al., 
2008; Simon et al., 2009; Xiang et al, 2009), the 
implications of glacier retreat on the microbial 
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communities in glacier-fed streams — most of them 
contained in benthic biofilms — remain elusive. 
Filling this knowledge gap is important, given that 
biofilms, even in glacier-fed streams (Battin et al., 
2004), orchestrate numerous ecosystem processes 
and contribute to large-scale carbon fluxes (Battin 
et al., 2008). Biofilms are remarkably diverse, partly 
because they assemble from different sources of 
biodiversity (for example, soil and groundwater) 
within the catchment (Hullar et al., 2006; Besemer 
et al., 2012). Species sorting from these sources can 
occur through local environmental conditions 
and is increasingly understood as a key mechanism 
underlying biofilm community assembly (Besemer 
et al, 2012). 

As glaciers recede, they change the hierarchical 
habitat template of the riverine landscape where 
environmental processes operating at local and 
regional scales differentially affect life in the 
glacier-fed streams (Ward et al., 2002; Brown et al., 
2007a; Milner et al., 2009). Melting glaciers mobilise 
ice-locked organic matter with implications for 
downstream carbon cycling (Singer et al., 2012) 
and heterotrophic activity (Milner et al., 2009). To 
explore how biofilm biodiversity and community 
composition in glacier-fed streams might respond to 
environmental changes as induced by glacier retreat, 
we surveyed 26 glaciers and their streams in the 
Austrian Alps. Substituting space for time, we 
investigated the relationship between microbial 
community composition and environmental vari- 
ables apparently linked to glacial retreat. Alluding 
to studies on invertebrate communities (Brown 
et al, 2007a; Milner et al, 2009; Finn et al., 2010), 
we hypothesise an increase of a diversity — the local 
richness in a given glacier-fed stream — with increas- 
ing deviation of the stream environment from a 
typical glacier-fed stream. Along the same line, we 
also postulate a decrease of p diversity — the varia- 
tion of communities among glacier- fed streams. 



Material and methods 

Study sites and sample collection 
The study sites covered 26 glaciers located along the 
main chain of the Austrian Alps (Silvretta, Otztaler 
Alps, Stubaier Alps, Zillertaler Alps, Venediger 
Group, Granatspitz Group, Glockner Group, 
Hochalmspitze and Ankogel; geographically ranging 
from 10.162-13. 278°E, 46.778-47.135°N and 
2100-2880 m above sea level; Supplementary 
Figure 1, Supplementary Table 1). Sampling was 
performed within 45 days in July and August 2010, a 
period with almost completely snow-free glaciers. We 
sampled glacier- fed streams within 10 m downstream 
of the glacier terminus to capture the physicochem- 
ical signatures from glacial and non-glacial sources 
upstream in the catchment (Ward et al., 2002; Brown 
et al., 2007a). Benthic biofilms were obtained by 
rigorously shaking ~ 80 stones (diameter ~ 2 cm) 



successively in sterile 50 ml tubes. This field proce- 
dure removed most of the biofilm from the stones. 
The supernatant was then filtered (0.2 urn filters, 
GSWP Millipore, Billerica, MA, USA) using sterile 
equipment. Two to three litres of streamwater were 
filtered through sterile 0.2-Lim filters (GSWP) and all 
filters were immediately frozen in liquid nitrogen in 
the field. Streamwater pH, electrical conductivity 
and water temperature were measured in situ using 
WTW probes (pH320, Cond340i). 

As a putatively prominent source of microbial 
diversity to the stream, we also sampled subsurface 
ice (0.3— 0.5 m depth) from 14 randomly chosen sites 
across the ablation zone of each of the 26 glaciers. 
We first removed the upper ice layer (~ 20— 30 cm) 
using sterile (ethanol-flamed) ice picks and then 
collected ice from beneath into sterile Whirl-Paks 
(Nasco, Salida, CA, USA). These samples were 
transported frozen to the next base camp. There, 
ice was melted within 1 hour and ~10l of icemelt 
were filtered through sterile 0.2-Lim filters (GSWP), 
which were immediately frozen in liquid nitrogen. 



Streamwater dissolved organic carbon and nutrients 
Sreamwater was analysed for dissolved organic 
carbon (DOC) concentration using a Sievers 900 
TOC Analyser (GE Analytical Instruments, Boulder, 
CO, USA) operated with an inorganic carbon removal 
unit. Prior to injection, DOC samples (GFF-filtered, 
pre-combusted) were automatically acidified in the 
analyser as recommended by the manufacturer. We 
determined a method detection limit of the Sievers 
900 according to US EPA guidelines and found it 
typically <6ugCl _1 . The detection limit was calcu- 
lated as three (approximate critical z-value for 99.9% 
confidence)-times the standard deviation of nine lab 
replicates of low DOC water (average DOC concentra- 
tion: 29.71 ± 1.43 ngCl 1 ). Concentrations of N-NH 4 , 
N-N0 2 and N-N0 3 in the streamwater were deter- 
mined using Continuous Flow Analysis (Alliance 
instruments, Salzburg, Austria). For further analyses, 
we used 'total nitrogen', which represents the 
summation of N-NH 4 , N-NQ 2 and N-NQ 3 . 



DNA extraction and PCR amplification 
Genomic DNA was extracted using the PowerSoil 
DNA extraction kit (MoBio, Carlsbad, CA, USA) 
according to the manufacturer's recommendations. 
The V4-V5 regions of the 16S rRNA gene was 
amplified in 25-ul PCR reactions containing 
0.5|imoll _1 of each primer (Thermo Scientific, 
Waltham, MA, USA), O^mmoll-MNTPs (Thermo 
Scientific), 40 ug bovine serum albumin (Thermo 
Scientific), 4mmoll _1 MgCl 2 (Thermo Scientific), 
1 U Taq-DNA Polymerase with the recommended 
PCR buffer (Thermo Scientific) and 4 ul DNA extract 
(2— 4ng DNA). Primers used for amplification of the 
16S rRNA gene were the universal 515F 5'- 
GTGNCAGCMGCCGCGGTAA-3' and 926R 5'-CCGY 
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CAATTYMTTTRAGTTT-3' (Quince et al., 2011), 
containing the 454 Titanium A and B adaptors, 
respectively. For each sample, two different sample- 
specific barcodes contained in the forward primer 
were employed to reduce barcode-specific bias 
(Berry et al, 2011). 

Samples were amplified using an initial denatur- 
ing step of 2min at 94 °C, followed by 30 cycles of 
30 s denaturation at 94 °C, 30 s annealing at 56 °C, 
1 min elongation at 72 °C and a final elongation for 
lOmin at 72 °C. Each PCR reaction included a 
negative control. PCR products were run on a 1% 
agarose gel and purified using the Gel Extraction Kit 
(Qiagen, Hilden, Germany). The purified PCR 
products were quantified by gel using the Gel Doc 
XR+ System (BioRad, Hercules, CA, USA) and 
pooled equimolar for pyrosequencing. 



454-pyrosequencing 

Amplicons were sequenced on a GS FLX Titanium 
Sequencer in Liverpool (Centre for Genomic 
Research, University of Liverpool, UK). Raw output 
files were filtered and de-noised using the software 
package AmpliconNoiseVl.O (Quince et al., 2011). 
After pre-clustering the sequences with PyroNoise 
(AmpliconNoiseVl.O), PCR single-base errors were 
corrected via SeqNoise (AmpliconNoiseVl.O), a 
sequence-based clustering method that performs 
alignment of the sequences. Chimeras were finally 
identified and removed with the Perseus algorithm 
at an intercept of a = — 7.5 and a coefficient of ft= 0.5 
(Quince et al., 2011). This filtering and de-noising 
steps reduced the number of 658 592 flowgrams to 
384 320 reads (read length 400 bp). A complete 
linkage algorithm on a 97% identity level clustered 
these reads to operational taxonomic units (OTUs) 
(AmpliconNoiseVl.O) and taxonomic assignments 
were determined by using a naive Bayesian rDNA 
classifier (Ribosomal Database Project; Wang et al., 
2007) at a 70% confidence threshold. 



Data analysis 

Dissimilarity matrices of the community composition 
were calculated by applying the relative abundance- 
based Bray— Curtis dissimilarity index. Based on 
these matrices, non-metric multidimensional scaling 
analysis was used to visualise differences in com- 
munity composition, and a non-parametric permuta- 
tional analysis of variance Anderson (2001) to test 
differences in community composition among habi- 
tats (using the R-packages 'vegan' and 'ellipse', 
Oksanen et al., 2011; Murdoch and Chow, 2007). A 
concomitant test for dispersion homogeneity (the 
multivariate analogue of variance homogeneity) is 
provided by the analysis of p diversity (see below). 

A principal component analysis was performed 
to illustrate the relationship and importance of 
z-transformed environmental variables. Further, 
after Hellinger-transformation of the microbial 
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community composition data, we performed a 
redundancy analysis to determine a potential effect 
of environmental variables on community composi- 
tion (Legendre et al., 2011; using the R-package 
'vegan': Oksanen et al., 2011). The explanatory 
values of environmental variables were further 
explored by running a forward selection test 
included in the 'packfor' package (Blanchet et al., 
2008). Additionally, these variables were tested for 
correlations with the relative abundances of the 
most abundant bacterial families. These Pearson 
correlation coefficients were visualised in a 
heat map. 

We used Diversity Estimation software (Quince 
et al., 2008), to estimate the 'true richness' of the 
communities by fitting the 454-pyrosequencing- 
obtained OTU abundances to a Sichel distribution 
(Sichel, 1974). For comparison, all samples were 
rarefied to the lowest number of reads obtained from 
an individual sample (1007 reads; vegan: Oksanen 
et al., 2011). Using multiple linear regression, we 
tested for correlations between environmental vari- 
ables and microbial richness. 

To estimate the probability that a biofilm commu- 
nity represented a random sample of the respective 
streamwater community, we performed a random 
resampling procedure on the samples from each 
glacier, using functions of the R-packages vegan, 
ecodist and gdata (Goslee and Urban, 2007; Oksanen 
et al., 2011; Warnes et al, 2011; Besemer et al., 
2012). Each tested sample pair consisted of the 
biofilm and the suspended community from one 
glacier-fed stream. Individual reads were sampled 
from the suspended community (the probability of 
each OTU to be sampled being its relative abun- 
dance) with replacement until the number of OTUs 
in this randomly assembled community equalled 
the richness of the respective biofilm community. 
This procedure was repeated to yield 1000 random 
assemblages. The probability of the biofilm commu- 
nity to fall within the distribution of these random 
assemblages was calculated as the percentage of the 
distances of the random assemblages to their 
centroid, which were as high as or higher than the 
distance of the biofilm community to the centroid. 
As a conservative approach, the biofilm community 
data set was reduced to those OTUs, which also 
occurred in the respective streamwater community. 

Dispersion is an expression of average community 
dissimilarity, equivalent to average Bray— Curtis 
dissimilarity, and serves as a measure of spatial 
community variation or p diversity (Anderson et al., 
2006a). Dispersion was computed as the average 
distance of communities to the respective group 
(that is, the habitat) centroid in the ordination space 
resulting from a principal coordinate analysis (aka 
metric scaling) of the Bray-Curtis dissimilarity 
matrix (Anderson, 2006b). Differences in dispersion 
among habitats were tested using a permutation 
procedure (Anderson, 2006b) with functions of 
vegan in R (Oksanen et al., 2011). 



The ISME Journal 



Biodiversity in glacier-fed streams 

L Wilhelm ef al 



1654 



The relationship between selected environmental 
variables (for example, streamwater temperature) 
and p diversity of streamwater and biofilm commu- 
nities was tested by multiple linear regression on 
distance matrices (Lichstein, 2006; Ptacnik ef al., 
2010). For this, the Bray— Curtis dissimilarity matrix 
was unfolded to a vector of pairwise dissimilarities 
between any two communities and analysed as 
dependent on pairwise average temperature and 
pairwise temperature difference in a multiple linear 
regression model. In this analysis, the hypothesis of 
interest concerns the effect of average temperature 
on dissimilarity, but temperature differences must 
be included as a covariate to account for environ- 
mental differentiation simultaneously affecting p 
diversity; significances are computed by permuta- 
tion (see Supplementary Information for details). 
Additionally, as a graphical analysis, we adapted a 
local polynomial regression fitting procedure 
(LOESS; Cleveland ef al, 1992) to estimate a 'local' 
P diversity at a given temperature of interest. Based 
on a neighbourhood-approach, we identified subsets 
of sites with similar temperatures, whose commu- 
nities were then included in the computation of 
P diversity as an average dissimilarity in a weighted 
manner (see Supplementary Information for details). 
In addition to temperature, these analyses were 
performed on elevation, electrical conductivity and 
pH, which were previously found to be significantly 
correlated to streamwater and biofilm community 
compositions. 

Data analysis and visualisation was performed in 
SPSS (IBM, Armonk, NY, USA), SigmaPlot and R 
2.13.0 (R Development Core Team, 2011). 

Results 

Microbial community composition and taxonomy 
The 454-pyrosequencing data set consisted of 
3489 ±1493 (average ± s.d.) reads per sample, which 
clustered into 11032 different OTUs. The original 
pyrosequencing output files are available at the NCBI 
Sequence Read Archive under the accession number 
SRX196420. Of all found OTUs, 2.9, 53.6 and 13.2% 
were unique to glacial ice, streamwater and biofilms, 
respectively (Supplementary Figure 2). In total, 7.4% 
of all OTUs were found in all three habitats, 19.8% 
were shared by biofilm and streamwater commu- 
nities and 0.6% occurred exclusively in biofilms and 
glacial ice. Non-metric multidimensional scaling 
revealed differences in community composition 
among all three habitats (Figure 1), which were 
further identified as significant by permutational 
analysis of variance (pseudo-F= 11.326, d.f.! = 2, 
d.f. 2 =75, P< 0.001). However, these differences in 
community compositions were confounded with 
differences in dispersion among habitats (P diversity, 
see below). While glacial ice communities were 
relatively constrained with little variation among 
glaciers, communities in the streamwater and bio- 
films varied broadly and overlapped (Figure 1). 



$ 0 

Q 



-2 




■ Ice 

▼ Streamwater 
• Biofilm 



-2-10 1 2 

NMDS1 

Figure 1 Community composition in glacial ice, streamwater 
and biofilms. Non-metric multidimensional scaling analysis of 
microbial communities based on Bray-Curtis dissimilarities in 26 
glacial ecosystems (stress = 0.16). Ellipses represent the 95% 
confidence interval. 



Taxonomic analyses of our sequence data revealed 
that Proteobacteria, Bacteroidetes and Actinobacteria 
were the dominant phyla in ice, streamwater and 
biofilms, whereas Cyanobacferier/chloroplasts were 
abundant only in glacial ice and biofilm communities 
(Figure 2a). Several less abundant phyla as Verruco- 
microbia or Nitrospira were primarily found in the 
streamwater, while they were rarely or not detected in 
biofilms and ice, respectively. Comamonadaceae, 
Chitinophagaceae and Oxalobacteraceae were the 
most abundant families detected in all three habitats, 
while Sphingobacteriaceae and Acetobacteraceae 
were primarily identified in the ice (Figure 2b). 
Streamwater communities were characterised by 
numerous families, many of which occurred neither 
in ice nor in biofilms (Figure 2b, Supplementary 
Figure 2). 



Relationship between environmental variables and 
microbial community composition 
To identify the potential effect of streamwater 
geochemistry and glacial characteristics on micro- 
bial community composition and diversity in gla- 
cier-fed streams, we first tested the relationship of 
the variables among study sites. Principal compo- 
nent analysis revealed streamwater electrical con- 
ductivity and pH to form a major environmental 
gradient across the glacier-fed streams. This gradient 
was also associated with differences in elevation 
and glacial coverage (% glaciation of the catchment) 
of the respective glaciers, and all four variables 
together accounted for most of the environmental 
variation among glacier-fed streams (Supplementary 
Figure 3; Supplementary Table 2). 

Using redundancy analysis, we found that the 
physicochemical environment significantly explai- 
ned variation in community composition in 
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Phyla 



Biofilm & Ice 

Cyanobacteria & Chloroplasts 
Deinococcus-Thermus 



Shared 

Proteobacteria 
Bacteroidetes 
Actinobacteria 
Acidobacteria 
Firmicutes 
Gemmatimonadetes 
Planctomycetes 



Biofilm 




Stream water 

Verrucom icrobia 
Nitrospira 
Chloroflexi 
OD1 

Euryarchaeota 
TM7 

Crenarchaeota 

OP11 

WS3 

OP10 

Chlamydiae 

Fusobacteria 

Spirochaetes 



Streamwater 



50 most abundant families 



Shared 

Comamonadaceae 
Chitinophagaceae 
Oxalobacteraceae 
Burkholderiales incertae sedis 
Xanthomonadaceae 
Sphingomonadaceae 
Cytophagaceae 
Gemmatimonadaceae 
Planctomycetaceae 
Pseudomonadaceae 
Gp3 

Caulobacteraceae 
Cystobacteraceae 



Ice & Streamwater 

Clostridiaceae 
Porphyromonadacea 
Sinobacteraceae 




Ice 

Chloroplasts (Streptophyta 

Chlorophyta) 

Sphingobacteriaceae 

Acetobacteraceae 

Moraxellaceae 

Microbacteriaceae 

Intrasporangiaceae 

Gp1 

Deinococcaceae 

Bradyrhizobiaceae 

Nocardiaceae 



Biofilm 



Biofilm 

Chloroplasts (Bacillariophyta 
Nocardioidaceae 



Biofilm & Streamwater 

Flavobacteriaceae 

Methylophilaceae 

Hydrogenophilaceae 

Rhodobacteraceae 

Bacillaceae 



Streamwater 



Streamwater 

Rhodocyclaceae 

Holophagaceae 

Nitrosomonadaceae 

Bdellovibrionaceae 

Opitutaceae 

Neisseriaceae 

Subdivisions genera incertae sedis 

Nitrospiraceae 

OD1 genera incertae sedis 

Anaerolineaceae 

Bacteriovoracaceae 

Desulfobulbaceae 

Helicobacteraceae 

Gp6 

Ruminococcaceae 
Legionellaceae 



Figure 2 Distribution of taxonomic groups in glacial ice, streamwater and biofilm. The percentage of (a) phyla and (b) the 50 most 
abundant families associated with each habitat is visualised in ternary plots. The position in the triangle indicates the relative abundance 
of each taxon among the three habitats; the size of the circle represents the relative abundance of taxa. 



streamwater [F= 1.399, P<0.01, n = 26) and in 
biofilms (F= 1.419, P<0.01, n = 2B). In the next 
step, forward variable selection found that biofilm 
community composition was significantly 
explained by streamwater electrical conductivity 
(P<0.05), pH (P<0.05) and streamwater tempera- 
ture (P<0.05). Likewise, electrical conductivity 
(P<0.01) and pH (P<0.01) were related to stream- 
water community composition, whereas stream- 
water temperature was not. Elevation explained 



variation in community composition of streamwater 
(P<0.01), but not of biofilms. Glacial coverage was 
neither related to biofilm, nor to streamwater 
communities (data not shown). 

Variation in electrical conductivity and pH corre- 
lated with the relative abundance of microbial 
families and to some extent even of phyla in the 
biofilms (Figure 3). For instance, Actinobacteria, 
Nitrospira and Verrucomicrobia varied with elec- 
trical conductivity, whereas Acidobacteria, 
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Figure 3 Heat map of Pearson correlations between environmental variables and families detected in biofilms. Colours represent the 
r-values of Pearson correlations between relative abundances of the 45 most abundant bacterial families and environmental parameters. 
Families were ordered according to taxonomic affiliations. We found several families that were strongly correlated to pH and electrical 
conductivity. Positive correlations with DOC were only found with Proteobacteria. Streamwater temperature and total nitrogen were only 
weakly related to the abundance of taxonomic families. 



Gemmatimonadetes and Proteobacteria were related 
to changes in pH (Supplementary Table 3). 

To better understand these differing community 
composition patterns, we looked for indications of 
species sorting as a possible mechanism of commu- 
nity assembly. We tested whether biofilms were a 
product of purely stochastic immigration from the 
source community suspended in the streamwater, 
by comparing biofilm communities with random 
samples of the streamwater communities. Resam- 
pling revealed that in all 26 glacial ecosystems, the 
biofilm communities significantly differed from the 
centroid of the random assemblages (probability of 
the biofilm community to fall within the distribu- 
tion of the random assemblages P< 0.001; 
Supplementary Figure 4). 



Relationship between environmental variables and 
microbial biodiversity 

a diversity was highest in streamwater, intermediate 
in biofilms and lowest in glacial ice, as revealed by 
'true richness' (Quince et al., 2008) and rarefied 
richness (Figure 4a). Based on our results from 
community composition, we investigated the relation- 
ship between elevation, electrical conductivity, pH, 
streamwater temperature and a diversity. Multiple 
linear regressions revealed elevation and electrical 
conductivity as main variables explaining the 'true 
richness' in streamwater (adjusted 1^ = 0.67, 



F= 10.508, P<0.001, 22 = 26) and biofilms (adjusted 
1^ = 0.43, F= 3.883, P<0.05, n = 26). Using simple 
linear regression, it was found that richness of both 
streamwater and biofilm communities decreased 
with higher elevation, though for biofilms this 
relationship was weak (Figures 4b and c); electrical 
conductivity was not significantly related to the 
richness of streamwater communities, while biofilm 
richness decreased with electrical conductivity 
(r 2 = 0.18, P<0.05, n = 26; data not shown). The 
simultaneous decrease of biofilm biodiversity with 
elevation and electrical conductivity seems counter- 
intuitive, as a higher contribution of glacial melt- 
water to streamflow results in lower electrical 
conductivity. This pattern is likely driven by an 
imprint of geology, particularly the prevalence of 
carbonate minerals, among our study streams, as 
indicated by a significant correlation between 
electrical conductivity and pH (r=0.56, P<0.01, 
n=26). 

In agreement with patterns of community compo- 
sition, (3 diversity among glacial ecosystems was 
highest in biofilms, intermediate in streamwater and 
lowest in ice (Figure 5a). We found a significant 
decrease of biofilm p diversity with streamwater 
temperature (Figure 5b); electrical conductivity, pH 
and elevation were not significantly correlated with 
biofilm p diversity (Supplementary Figure 5). Micro- 
bial p diversity in streamwater did not correlate with 
any of the tested variables (data not shown). 
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Figure 4 Microbial a diversity, (a) Boxplot representation of 
rarefied richness (all samples rarefied to 1007 reads) and 'true 
richness' estimated by fitting Sichel distribution curves to the 
abundance distributions obtained from the 454-pyrosequencing 
data. Median (line), 1st and 3rd quartile (box margins), 5 and 95% 
percentiles (whiskers), outliers (points). Different letters on each 
box represent significant differences between habitats as deter- 
mined by Wilcoxon tests followed by Bonferroni correction 
(P<0.01). (b,c) Microbial richness decreased with higher eleva- 
tion. Linear regression include 95% confidence intervals in the 
(b) streamwater and (c) biofilms. 



Discussion 

Glaciers retreat worldwide, making it necessary to 
better understand the impacts of retreating glaciers 
on glacial ecosystems. Our study streams 
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Figure 5 Spatial variation of microbial communities, (a) f) 
diversities of the three habitats as the average distance of 
communities to the group (that is, habitat) centroid after principal 
coordinate analysis of the Bray-Curtis dissimilarity matrix. 
fS diversity differed significantly among habitats (pseudo- 
F=17.082, d.f. a = 2, d.f. 2 = 75, P<0.001 for all pairwise compar- 
isons), (b) Graphical analysis of biofilm [5 diversity as a function of 
streamwater temperature by local polynomial regression fitting. 
Multiple linear regression on distance matrices revealed p diversity 
to decline significantly with increasing streamwater temperature 
(pseudo-f= —4.77, P<0.05). The shaded area gives a bootstrap 
confidence interval for the generated trendline. See Materials and 
methods and Supplementary Information for computational details. 



encompassed a gradient of elevation, glacial cover- 
age and streamwater electrical conductivity, suggest- 
ing distinct hydrological settings of the glacier-fed 
streams at the catchment scale as a result of variable 
contributions of icemelt, snowmelt and groundwater 
to streamflow (Brown et ah, 2003; Brown et ah, 
2007a; Milner et ah, 2009). Additionally, catchment 
geology and hydrochemical processes may imprint 
the observed environmental gradient, as indicated 
by the positive correlation of streamwater electrical 
conductivity and pH. At local scale, depending on 
geomorphology (for example, slope and slope 
breaks), the ablation zone of some of our study 
glaciers was fragmented following thinning, and 
meltwater therefore flowed over bare rock exposed 
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to solar radiation (Supplementary Figure 6). This 
caused elevated and outlying streamwater tempera- 
tures in these systems and weakened the relation- 
ship between streamwater temperature and glacial 
coverage across all sites. The physicochemical 
signatures in glacier-fed streams thus reflected 
catchment scale and local processes operating along 
the various flow paths (that is, supraglacial, engla- 
cial, subglacial and non-glacial). 

The most abundant phyla identified in our study 
streams are typical for freshwater ecosystems 
(Tamames et al., 2010). Their regular occurrence in 
other ice ecosystems (Simon et al., 2009; Xiang et al., 
2009; Anesio and Laybourn-Parry, 2012) suggests that 
they comprise cryophilic taxa, as also supported by 
cultivation-dependent approaches (Cheng and Foght, 
2007; Loveland-Curtze et al., 2009). Streamwater 
electrical conductivity and pH were associated with 
community composition in biofilms at the level of 
OTUs, family and even of phyla (Figure 3, 
Supplementary Table 3). These results are in line 
with soil studies (Fierer and Jackson, 2006) suggest- 
ing that pH — a prime physiological control on single- 
celled organisms — structures microbial communities 
even at the phylum level (for example, Acidobacteria, 
Proteobacteria), and with studies on invertebrate 
communities in glacier-fed streams (Brown et al., 
2007b). Certain biofilm taxa, some of them (for 
example, Nitrospira) having critical roles in biogeo- 
chemistry, seem particularly prone to physicochem- 
ical shifts. Streamwater temperature, known to 
influence benthic microbial community composition 
in streams (Hullar et al., 2006), appeared to further 
affect biofilm community composition, probably at a 
more local scale given the pattern of streamwater 
temperature among our study streams. Expected long- 
term changes of physiochemical conditions following 
glacial retreat could potentially have widespread 
implications for glacial stream ecosystem functions. 
Owing to their sessile mode of life, biofilm micro- 
organisms are likely even more susceptible to 
such changes compared with transient cells in the 
streamwater. 

Despite the cold and mainly oligotrophic condi- 
tions in glacial ecosystems, our 'true richness' 
estimates for streamwater and biofilm communities 
were comparable to those from non-glacial streams 
(Besemer et al., 2012). In line with the present study, 
Besemer et al. (2012) reported higher microbial 
richness in streamwater than in biofilm commu- 
nities. We attribute the differences in a diversity 
between ice, streamwater and biofilms to various 
degrees of environmental harshness [sensu Jacobsen 
and Dangles, 2011) and metacommunity (that is, a 
set of local communities linked by dispersal) size 
(Leibold et al., 2004), which interact to determine 
biodiversity. In fact, glacial ice is a harsh but 
comparatively constant environment (Hodson 
et al., 2008) harbouring a constrained microbial 
community, which may be particularly well adapted 
to its environment; cell immigration is assumedly 



dominated by atmospheric deposition (Hervas and 
Casamayor, 2009). Glacier-fed streams, however, are 
dynamic with pronounced temporal fluctuations 
(Brown et al, 2007a; Milner et al., 2009), and collect 
microorganisms from various glacial (for example, 
subglacial, englacial and supraglacial runoff; Anesio 
and Laybourn-Parry, 2012) and non-glacial sources 
(for example, groundwater; Brown et al., 2007a; 
Milner et al., 2009). In fact, microorganisms from 
cryoconite holes (that is, holes with high microbial 
activity at the glacier surface, which form around 
solar-heated debris; Edwards et al., 2011), the 
englacial environment (Anesio and Laybourn- 
Parry, 2012), and also from groundwater, atmo- 
spheric deposition (Hervas and Casamayor, 2009), 
residual snow (Hervas and Casamayor, 2009) and 
adjacent soils and rocks (Bardgett et al., 2007; 
Schiitte et al., 2010) may form a metacommunity 
that potentially contributes to the community in the 
glacier-fed streams. We suggest that decreasing oc 
diversity in the streamwater with elevation may be 
due, at least in part, to less diverse sources of 
microorganisms upstream in the catchment. Appar- 
ently glacial ice contributed only marginally to the 
microorganisms in the streamwater, whose a diver- 
sity may therefore increase as other compartments of 
the metacommunity (for example, groundwater, 
soils) gain importance. Similarly, the low common- 
ality of taxa both at the OTU (Supplementary 
Figure 2) and family level (Figure 2b) between ice 
and biofilm communities suggests minor contribu- 
tions of the ice communities to biofilm assembly in 
the glacier- fed streams. 

Differing patterns of community composition and 
diversity suggest deviating assembly mechanisms in 
biofilm and streamwater communities. In the frame of 
metacommunity theory (Leibold et al., 2004), species 
sorting, where the local environment and biotic 
interactions select from the metacommunity, 
would be a candidate mechanism for biofilm comm- 
unity assembly (Besemer et al., 2012). In contrast, 
streamwater communities could likely be explained 
by mass effects, determined by large cell influx rates 
and short residence times limiting the influence of 
environmental factors. The local environment in the 
glacier-fed streams would thus select microorganisms 
from the streamwater for biofilm formation, the 
diversity of which is influenced by the size of the 
metacommunity upstream. The latter may be affected, 
for instance, by the diversity of hydrological flow 
paths upstream and the various habitats they connect, 
and also by a possible catchment scale imprint on 
geochemistry (for example, electrical conductivity, 
pH). We found that stochastic immigration of OTUs 
from the streamwater into the biofilms was unlikely to 
explain the observed community composition of 
biofilms, thus supporting the assumption that species 
sorting has a role in biofilm community assembly. 
Collectively, these findings suggest that the various 
components of the hierarchical habitat template 
differentially influence the composition of microbial 
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communities in glacier-fed streams, which complies 
with general stream ecology (Ward et al, 2002). 

Ice, streamwater and biofilm communities exhib- 
ited different levels of spatial variation among all 
study sites, likely reflecting community responses to 
various degrees of environmental variation in these 
habitats, and different assembly mechanisms. The 
relatively small range of streamwater temperature 
was apparently sufficient to influence the spatial 
variation of biofilm communities among glacier-fed 
streams. We suggest therefore that sensitive warm- 
ing of glacier-fed streams, driven by streamflow 
shifts from glacial to non-glacial sources or local 
flow over exposed bedrock in a fragmenting 
glacial landscape, could have an impact on the 
distribution of microbial diversity in these systems 
(Vincent, 2010). 

Substituting space for time, our findings suggest 
biofilms in glacier-fed streams as sentinels of glacier 
retreat and expand the current knowledge on the 
threatened biodiversity in the European Alps (Brown 
et al, 2007a; Finn et al, 2010; Jacobsen et al, 2012). 
Invertebrate a diversity is thought to increase as 
glaciers recede, because of species migrating 
upstream (Brown et al., 2007a; Milner et al, 2009; 
Finn et al, 2010). Microorganisms, however, primar- 
ily disperse downstream with water flow. Therefore, 
we suggest that both shifts in the metacommunity 
upstream of the glacier-fed stream and changes of the 
local physicochemical environment, as induced by 
glacial retreat, alter biodiversity and composition of 
the streamwater communities and, via species sorting, 
also of biofilm communities. Glacier retreat may 
increase microbial a diversity in glacier-fed streams, 
while concomitantly reducing (3 diversity, and 
thereby contributing to the homogenisation of biofilm 
communities among glacier-fed streams — similar to 
the patterns observed for invertebrates (Jacobsen and 
Dangles, 2011; Jacobsen et al, 2012). Spatially 
isolated and sporadically occurring invertebrate 
populations face an elevated risk of local extinction 
as glaciers recede, which ultimately reduces their 
spatial variation (Finn et al, 2012; Jacobsen and 
Dangles, 2011; Jacobsen et al, 2012). The responsive- 
ness of biofilms to the local environment in glacier- 
fed streams invokes a similarly patchy distribution 
susceptible to homogenisation following glacier 
retreat. These findings call for more research explor- 
ing the implications of microbial biodiversity shifts in 
glacier-fed streams for ecosystem processes and 
downstream biogeochemistry. 
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